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The dynamic resistance of a sphere with a general inhomogeneous slip boundary condition is analysed 
in Newtonian unbounded uniform flow at low Reynolds number. The boundary condition is treated as 
a perturbation to a homogeneous sphere, assuming that the slip length magnitude b is much smaller 
than the sphere radius a. To first order, the effect of inhomogeneous slip is the same as that of a radial 
deformity of magnitude b. Full resistance tensors are presented and the dynamics of a hemispherical 
inhomogeneous sphere, such as a Janus particle, are explicitly calculated. 

The motion of micro- and nanospheres in fluids is of increasing importance in several emerging fields 

of research. Spheres have attracted interest for biomedical applications such as targeting cancer cells 

{l|, Q] . Magnetic microspheres have frequently been used in recent research [3:] , while small particles have 

the potential to play a significant role in microfluidic and MEMS technologies. Fabrication of micro- and 
> 

' nanospheres is widespread, because it is relatively easy to obtain spherical particles with homogeneous 
surface chemistry due to minimization of surface tension. Inhomogeneous spherical particles such as 
dual- hemispherical Janus particles [J] (including magnetic Janus particles |^) are also being fabricated. 

Manipulation of small particles in fluids has been achieved by various means, particularly using 
electromagnetic methods, which have been described analytically [gI. Other factors relevant to particle 
manipulation that are less well studied include orientation dependence |l|, particle inhomogeneity and 



1 



surface slip. Slip occurs when there is non-zero relative motion between a solid surface and the fluid 
immediately adjacent to that surface, violating the non-slip boundary condition (NSBC). The presence 
of slip on a particle surface alters the shear forces on that surface, which are otherwise very stable at low 
Reynolds number, thereby introducing novel manipulation properties that do not require application 
of external fields. The rotational dynamics are likely to be of intrinsic importance for inhomogeneous 
and Janus particles. In the last 15 years, there has been great interest in practically achieving surface 
slip in microfluidic flows [7, 8]. Experimental work has shown that, for a Newtonian fluid at a smooth. 



hydrophobic surface, non-zero 



on the microscale 
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slip does occur and can be greatly intensified using surfaces engineered 



13l |. Complementary theoretical work has most frequently used Navier's 



formulation of the slip boundary condition 

Analysis of spherical and near-spherical particles in fluid flow at low Reynolds numbers has consis- 
tently been of interest over the past two hundred years. Happel and Brenner \l\ produced a summary 
of relevant analyses, covering the flows surrounding particles and the resulting dynamics. In particular, 
they presented the solution for a sphere with a homogeneous slip boundary condition. The problem of 
a slightly deformed sphere with a homogeneous slip boundary condition has since been solved approxi- 
mately 

In this Letter, the dynamic response of an impermeable spherical particle with a spatially varying, 
inhomogeneous slip boundary condition is analysed for the first time in streaming Newtonian flow at low 
Reynolds number. Slip is formulated as a general function of surface area and applied as a perturbation 
to the solution for a homogeneous NSBC sphere. This approach gives first-order results for the resistance 
tensor of the particle for any slip boundary condition and any relative flow direction, providing the slip 
length magnitude is much smaller than the sphere's radius. Results are calculated for some specific 
inhomogeneous boundary conditions, including those relevant for Janus particles. 

The equation of motion for steady flow of an incompressible fluid (V.v = 0) at low Reynolds number 
is Stokes' equation, 

ryVV = Vp, (1) 
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where v is the velocity field for a Newtonian fluid of density p and viscosity 77, and the pressure term p 
absorbs any body forces. Solutions to Stokes' equation are dependent on flow geometry and boundary 
conditions. There is a general methodology for finding v for unbounded uniform fiow external to a 
particle in viscous fiow, and the force and torque on a spherical particle [l^. For inhomogeneous 
boundary conditions, it is not practical to analytically find the full solution, which requires several 
expansions of spherical harmonic functions. 

Navier's slip boundary condition \^ states that the component of fluid velocity tangent to a solid 
surface (wii) is proportional to the shear rate at the surface, 



;|| = 6n • ( Vv + (Vv)^ ) • (I - nn) . (2) 



Here, n is the unit normal to the surface and the equation is valid in the rest frame of the wall [Vi 
The parameter h has units of length and is often referred to as the 'slip length'. For an impermeable, 
homogeneous surface, Eq. [2] reduces to v\\ = where is the spatial gradient of v\\ normal to the 
solid surface and b\\ = 6 if the solid-liquid interface is assumed to be planar, as has been the case in most 
recent slip studies. The derivation of h for a curved surface yields 

111 

where R is the radius of curvature of the interface, defined as positive for a convex solid surface. Note 
that 6|| = when 6 = 0; 6|| k. h when b « R] and if R and h are of comparable magnitude, 6|| > h. 

In order to describe an inhomogeneous boundary condition on the surface of a sphere, the slip length 
hf (9, (p) is used, where 9 and (p are spherical polar coordinates and f {9,(f)) is a function of 0(1). When 
h is small compared with sphere radius a, a perturbed flow around the sphere can be expanded in powers 
of e = — (6||/a), where |e| « 1 and therefore b\\ ~ h: 

v = v(0)+ev«+eV2)+... (4) 
The inhomogeneous version of the flat-surface slip boundary condition (Eq. [2|) is 
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where / (6, 
conditions, 



V = -eaf {0,4>)^, 



(5) 



can be expanded in terms of surface spherical harmonics fk(0,(l)), such that f{6,(j)) = 
Comparing Eqs. S] and [5l the first two terms in the expansion of v give the boundary 



v(0) (r = a,e,, 



(7V 



(0) 



(6) 



dr 



For incident unbounded uniform flow characterised by the velocity vector U, the boundary conditions 
far from the sphere are 



(r = 00,^,0) = U, 
v(i) (r = oa,e,(t>) = 0. 

The solution for the zeroth-order fleld v*^") is the usual solution for a homogeneous sphere, 



(7) 



v(0) = U • 



3a /I rr 



4 \ V T'~ 



4 V r 



I rr 

3^ 



(8) 



The corresponding pressure, drag and torque are similarly well known. To find v^^^ , a spherical harmonic 
expansion is used. 



,(1) 



k=0 



(1) 
k ' 



(9) 



where each v^^-* {k = 0,1,2,...) satisfies Stokes' equations and (from Eqs. [U [7] and [8]) the boundary 
conditions 



v^j}^ (r = 00,6*,. 



(10) 



0. 
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The requirement that Vr {r = a,0,(j)) = was not exphcitly stated when determining the boundary 
conditions, but is met by Eqs. [TlandfTOl 

For a fuh description of the spherical harmonic formulation of v(^) and the relating dynamics, the 
reader is referred to the analogous treatment of a slightly deformed sphere given by Happel and Brenner 
[l^ . in which the sphere's surface takes the locus of points 

r = a{l + ef{e, (/))). (11) 

When this surface is treated as a perturbation to a sphere under the NSBC, the formulation is coincident 
with the first-order problem defined by Eqs. [H [H] and [71 with e being replaced by e. To first order in e, 
the force and torque on the inhomogeneous slip sphere are 



F = f(0) + eF(i) + O (e^) 

= 67r7?aU - Mnva ( U/o - - (U • V) V {r^f2)] + O (e^) , 



T = t(0) + eT(i) + O (e^) 

= + -levTT/a^U X V (r/i) + O (e^) . 

The dynamics for a rotating, slightly deformed sphere can be calculated following a similar method, 
once the boundary conditions are carefully checked. Again the perturbation of the velocity field (Eq. [3]) 
is used, where |e| « 1. The inhomogeneous slip boundary condition is 

dv 

V = cj X r-ea/(6',(^)— . (14) 
Comparing Eqs. H] and [T^ the boundary conditions at the sphere's surface are 



^(0) (^r = a,9,(j)) = u; X r 

{r = a,e,cf>) = -afiOA) 



(15) 
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The boundary conditions at r = oo are v^*^) = v^^^ = 0. These boundary conditions are also analogous 
to the spherical harmonic treatment of a slightly deformed sphere [ij]. Conversely to the case for 
translational motion, only odd harmonics generate a force on the rotating particle, while only even 
harmonics generate torque. To first order in e, 

F = + eTTTya^w X V (r/i) + O (e^) , (16) 



T = -87r7?a=^^ + (^^^ 24^r?a3 (ufo - 1 (c^ . V) V (r2/2)) + O (e^) . 



(17) 



Stokes' equations and the boundary conditions are linear, so solutions for translation and rotation 



can be superposed to give the full resistance tensors 
The first-order tensors about the centre of an inhomogeneous, slipping sphere are 



(18) 



'^if = ^^^/o - 3^^^^^- , (19) 



Cip =D'.f = -e.,,k'Vk{rfi), (20) 

f^S^ = -3 (S.^fo - ^ViVj (rV2)) , (21) 

where Sij is the Kronecker delta, £ijk is the Levi-Civita permutation symbol, and tensor components are 
given in a Cartesian frame. 

As noted above, the flow field and resulting dynamics for a slip length b are dynamically equivalent 
to surface deformation —ea for a sphere in viscous flow when e « 1 and e « 1. The mathematical 
equivalence of these approaches does not mean that they describe physically similar phenomena. The 
boundary conditions are coincident to flrst order, but not to second order, because the slip boundary 
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conditions for v^^^ and higher-order terms do not contain second and higher derivatives of the flow 
profile. When higher-order terms are considered, the difference between b and fty becomes important. 
Another physical difference is that surface deformation is independent of strain rate, whereas the rate 
dependence of b (and, therefore, the linearity of Eq. [2]) is an unresolved topic of active research ^j. 
For an explicit calculation, consider the specific inhomogeneous boundary conditions given by 

[7 {0<9 <e') 
f{e,cl>) = \ (22) 
[1 {6' <6 <tt). 

which are shown schematically in Fig. [H along with the Cartesian (x, y, z) and spherical polar (r, 9, (p) 
co-ordinate axes. To describe this boundary condition, an expansion of normalized Legendre polynomials 
Pn (cos 9) is used, 



f{0A)=f0 + fi+f2 + :. 

= - (1 cos 61' 7 (1 - cos 6*')) Po (cos 61) 

3 (23) 
+ - (cos^ 9' - 1) (1 - 7) Pi (cos 9) 

+ ^ (cos^ 9' - cos 9') (1 - 7) P2 (cos 9) + .... 
For the hemispherical boundary condition relevant to Janus particles {9' = 90°), using Eqs. 1121 1131 and 



F 



/ b»(l + -f)\ 
evrryaU (l - " , (24) 



o^ll /3r \ 

T = -evrr/a^-l (1 -7)U X V f — cos^J . (25) 

The rotational motion reaches equilibrium when the vector V(r/i) is parallel to U. In Cartesian 
coordinates, the torque about the centre of the sphere reflects the symmetry of the boundary condition 
about the z axis: 
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T = -^vTT^a^^ (1 - 7) {Uyit - C/.y) ■ (26) 

The dynamic response of a slipping sphere provides a mechanism for manipulation of that sphere 
when there is fluid motion relative to the particle. The sphere is lubricated, while asymmetry of the 
boundary condition relative to the flow direction allows manipulation of orientation. Such dynamics are 
important, for example, when the particle is tethered stationary relative to flow (but free to rotate), 
or competing with other applied forces. Application of external fields is not required. The effect 
of slip increases with (b/a), and it is important to recognise that the slip length is not in principle 
restricted by particle size. Practical manipulation of spheres using asymmetric slip will be dependent 
on controlling the shape of the particle and the method of inducing slip. Superhydrophobic surfaces 
generally have physical structure of 0(1) /im while providing Newtonian slip lengths up to 0(10) /xm, 
while slip lengths for water on smooth, hydrophobic surfaces are of 0(10) nm. Development of reliable 
experimental methods for slip length measurement is an area of active research ^. The major methods 
currently used do not incorporate significant surface curvature, nor do they precisely test the linearity 
of Navier's slip length. The calculations described here provide a method for measuring b. 

The influence of thermal motion on orientation effects can be estimated by considering the Peclet 
number for the rotating sphere. The mechanical energy scale is the work done by flow parallel to the x 
axis (Ux) when rotating an inhomogeneous slip boundary condition (7 = 1) sphere 90° about the y axis, 
so that the Peclet number is 

^ 9W (27) 

where T is temperature and ks is Boltzmann's constant. Figure O is an indicative plot of the conditions 
at which thermal motion becomes signiflcant (Pe ~ 1) for water at room temperature. 

To summarize, first-order solutions for low-Reynolds number sphere dynamics that are general with 
respect to boundary condition inhomogeneity and the relative orientation of unbounded uniform fiow 
have been derived. The full resistance tensors, inclusive of rotational motion, have been presented. 
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The slip perturbation is mathematically equivalent to a slightly deformed sphere, but the governing 
physical processes for these situations are significantly different. The calculations could be applied to 
novel applications utilizing inhomogeneous slip. They are just as likely to be useful for characterising the 
intrinsic properties of particles that are already in use. The analysis will be useful for those interested 
in slip length measurement, especially for smooth surfaces. Thermal fluctuations will play a significant 
limiting role in the usefulness of slip-based orientation techniques, especially in atmospheric conditions 
on the nanoscale. 
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FIG. 1. Figure 1(a) defines Cartesian (x,y,z) and spherical polar (r, 6',(/)) coordinate axes. The centre 



of the solid sphere of radius a considered in the analyses is located at the origin. Figure 1(b) shows a 
specific inhomogeneous slip configuration (Eg. I22p . 
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U = 1 mm s" 

U^ = ^ ixms'^ 

b|i;a = 0.1 

b„/a = 0.001 




1E-8 1E-7 1E-6 1E-5 1E-4 1E-3 
a / m 



FIG. 2. Plots of the locus of points at which Pe = 1 (Eg. 1270 . For each line, thermal motion is significant 
at points closer to the x axis. The fluid is assumed to be water at 300 K. 
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